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Analytical calculations and recent numerical experiments have shown that a sizable of the mass 
in our Galaxy is in a form of clumpy, virialized substructures that, according to jjl . can be as light 
as 1O~ 6 M0. In this work we estimate the gamma-rays flux expected from dark matter annihilation 
occurring within these minihalos, under the hypothesis that the bulk of dark matter is composed by 
neutralinos. We generate mock sky maps showing the angular distribution of the expected gamma- 
ray signal. We compare them with the sensitivities of satellite-borne experiments such as GLAST 
and find that a possible detection of minihalos is indeed very challenging. 
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Analyses of the anisotropics in the cosmic microwave 
background radiation |2| find that the matter density 
content of the universe is approximately six times larger 
than the baryonic one, in agreement with the observed 
abundance of light elements Q and the matter power 
spectrum inferred from galaxy redshift surveys Q . These 
observations involve very different physics, but unam- 
biguously indicate that the Universe contains a signifi- 
cant amount of non-baryonic cold dark matter (CDM). 
The nature of CDM is presently unknown. However, 
weakly interacting massive particles (WIMPs) are re- 
garded as generic particle candidates for CDM, as they 
naturally arise in extensions of the standard model of par- 
ticle physics. WIMPs are a particularly attractive CDM 
candidate, since the average mass density of a stable relic 
from the electroweak scale is expected to be close to crit- 
ical 5] . A further exciting feature of WIMPs is that they 
can be investigated in ongoing and future astrophysical 
direct @ and indirect searches @ and in upcoming lab- 
oratory experiments. 

Indirect CDM searches focus on measuring the dif- 
fuse flux of CDM annihilation products $ 7 = $ SUSY x 
(jjcosmo |gj rp ne p ar ticle physics dependence embedded 
in $ SUSY involves the WIMP's annihilation cross-section 
and branching ratios as well as the final photon en- 
ergy spectrum and can be calculated from the underlying 
WIMP field theory. In this work we assume the nearly 
mass independent best value for the annihilation cross 
section of 2 x 10 _26 cm 3 s _1 , which represents an upper 
bound still allowing the neutralinos to be the dominant 
CDM component. The geometry dependence of the flux 
is contained in $ cosmo that is found by integrating the 
square of the neutralino mass density along the line-of- 
sight and therefore is very sensitive to the presence of 
virialized structures representing maxima of the DM den- 
sity field. In this work we focus on the annihilation signal 



produced within our Galaxy under the hypothesis, sup- 
ported by several high resolution numerical experiments, 
that part of its mass is in the form of subhalos. 

CDM models are characterized by their excess power 
on small scales that leads to a logarithmic divergence of 
the linear density contrast at large wavenumbers A = 
Sp/p cx ln(fc). Recent analytical calculations have proved 
that A shows exponential damping for k > k cut = 
C(l)/pc. The cut-off scale fc cut is given by viscous (colli- 
sional) processes before and free (collisionless) streaming 
after kinetic decoupling at Tkd = O(10) MeV, both lead- 
ing to exponential damping of the linear CDM density 
contrast [3, Hoi ITlj . During the linear regime most of 
the statistical power typically goes in density contrasts 
A with k ~ fc cut . So fc cu t sets also the typical scale for the 
first halos. that typically form at a redshift z n = 60 ± 10 
(for the best fit WMAP matter density) [13,113, when 
the mass variance a = 1 on a comoving length scale 
R = C(l)pc. 

The fate of these early minihalos in the non linear 
regime when they typically merge into larger halos, rep- 
resenting higher levels in the hierarchy, can only be 
followed through numerical experiments. The limited 
dynamical range explored with currently available N- 
body codes makes it very challenging to perform nu- 
merical experiments with a resolution as high as k~^ t on 
a computational box large enough to guarantee statis- 
tical significance. Yet, [l| have recently simulated hi- 
erarchical clustering in CDM cosmologies up to a dy- 
namical resolution fc ros ~ fc cut in a small high resolution 
patch which is nested within a hierarchy of larger low 
resolution regions. They find a steep halo mass func- 
tion drc(M)/dln(M) oc M' 1 exp[-(M/M cut )~i], with 
M cut = 5.7 x 10 _6 /i _1 M Q cx fc c ~ t 3 in the mass range 
[10~ 6 , 10~ 4 lM fT) w hose extrapolation fits well that of 
Galactic halos ITJ. 
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Since |l| stopped their experiment at z = 26 and 
probed a small volume with a limited dynamical range, 
no general consensus exists on whether these early struc- 
tures survived within massive galactic halos until today 
and what their mass function and distribution is. In- 
deed, although their mass concentration is probably large 
enough for them to survive the gravitational tides expe- 
rienced in the external region of the massive host halo 
(see, however 0), they could be torn into tidal stream 
by close encounters with individual stars in the galaxy 
UJ. These issues could only be settled by simulating 
the gravitational clustering out to z = on a volume 
large enough to overlap the results with those of other, 
high resolution numerical experiments of Galaxy size ha- 
los f which is beyond current numerical limita- 
tions. Yet, [J, have suggested that these minihalos 
might significantly contribute to the total annihilation 
flux in our Galaxy. 

In this work we compute the expected 7-ray flux pro- 
duced by a population of sub-Galactic halos under the 
optimistic hypotheses of |lj that all minihalos of masses 
> lO _6 Af that survived gravitational tides populate the 
Galactic halo at z = 0, trace its mass, share the same 
self-similar cuspy density profile and have a steep mass 
function cx AT -1 up to > 1O 1O M0. It is worth stressing 
that our predictions represent an upper limit to the ex- 
pected photon flux and to its experimental detectability 
that we will investigate for the case of a satellite-borne 
experiment like GLAST 0. 

Our analysis extends those of [H E |H HJ H3 since 
we consider subhalos much lighter than 1O 6 M and since 
we use Monte Carlo techniques to explicitly account for 
the contribution of very nearby minihalos to the total an- 
nihilation flux. The relevant properties of our minihalo 
population are: a steep mass function dn(A/)/dln(Af) oc 
A/ -1 in the range [10 -6 , 10 10 ]Mq, a spatial distribution 
that trace the smooth mass com pon ent in our Galaxy 
and a NFW mass density profile [23] for both the MW 
and for all the subhalos: p x = p s (r/r ;! )~ 1 (l + r/r s )~ 2 , 
where both the scale radius r s and the scale density p s 
depend on the concentration parameter c which is a func- 
tion of mass, redshift and the cosmological model. Since 
we aim at exploring the most optimistic scenario for neu- 
tralino annihilation we do not consider here the results of 
some recent high resolution N-body experiments l2r| 
which suggest that halo density profiles near the center 
have a continuously varying logarithmic slope, and thus 
are significantly shallower than the cuspy NFW profile. 
In this work we assume the flat, ACDM "concordance" 
model and use the numerical routines provided by |24| to 
compute c. The smallest halos of mass 1O _6 M turn out 
to have c ~ O(40) at z = 0. 

The previous assumption allows one to specify the 
number density of subhalos per unit mass at a distance 
r from the Galactic Center (GC): 

C\/f \ Al\r-2 6( r — r min{M)) _i 3 

Psh (M, r) = AM ^ r/rMW){1 + r/rMW)2 M Q kpc 

(1) 



We set A such that 10% of the MW mass (A/mw — 
1O 12 M ) is distributed in subhalos with masses greater 
than lO 7 Af to match the results of As a result 

about 53% of the dark mass within our galaxy is not 
smoothly distributed but contained within ~ 1.5 x 10 16 
subhalos with masses larger than lO _6 Af , correspond- 
ing to ~ 100 pc -3 halos in the solar neighborhood. 

The term in Eq. Stakes into account the effect of 
gravitational tides which, according to the Roche crite- 
rion, disrupt all halos within r m i„(M) from the GC [27| 
in an orbital period. At z — r m i n is an increasing 
function of the subhalo mass implying that no subhalos 
survive within r m i„(10 -6 Af ) ~ 200 pc. 

In this work we wish to predict the annihilation sig- 
nal expected in a 100° x 100° field of view (f.o.v.) with 
an angular resolution of A£l = 10~ 5 sr, matching those 
planned for the GLAST experiment, in both the GC and 
the anticenter (AC) directions. This requires to evaluate 
the line-of-sight integral: 

$™(^,0)=/ dSl' f p 2 x (r(X,4,'))d\(r^'), 

JAn(4>,e) Ji.o.s 

(2) 

where ip is the angle-of-view from the GC, 8 is 
the angular resolution of the detector, p x (r) is the 
mass density. The distance from the GC r is r = 

\J \ 2 + Rq 2 — 2XRq cos -0, A is the distance from the ob- 
server and Rq — 8.5 kpc is that of the Sun from the GC. 

We evaluate of Eq. [3 in two steps. First, we nu- 
merically integrate Eq. [2 in which p x (r) is the sum of 
NFW profiles corresponding to all subhalos distributed 
according to Eq. ^ along the l.o.s.. This gives the 
average subhalo contribution to the Galactic annihila- 
tion flux within 10 _5 sr along the direction (ip,0). We 
found that the average contribution to <j) cosmo ; in units of 
GeV 2 cm~ 6 kpcsr, of subhalos in the range lO _6 Af < 
M sh < 10 6 M & is 4 x 10~ 5 at the GC, then slowly de- 
creases to 10~ 5 at ip = 50° and to 4 x 10~ 6 at larger 
angles. For (j> > 25° it dominates over the Galactic 
foreground accounting for neutralino annihilation in the 
smooth Galactic halo of 4.7 x lO n M . 

To estimate the variance to the average flux we use 
Eq. [21 to generate several independent Monte Carlo real- 
izations of the closest and brightest minihalos. For each 
subhalo mass, we only consider objects which are close 
enough to guarantee $ cosmo > l[r 6 GeV 2 cm -6 kpc sr. If 
no such halo exists then we still Monte Carlo generate the 
100 closest objects in that mass range, since we expect 
the bigger fluxes to come from the closer halos. 

The contribution of nearby structures to the annihila- 
tion flux can be appreciated from Fig^ which shows the 
sky distribution of the annihilation signal in one of our 
Monte Carlo realization in the direction of the GC. 

To compute the total annihilation signal from all sub- 
Galactic halos we have followed the same procedure as 
[2fj| and have generated several Monte Carlo realizations 
of subhalos with 10 6 M Q < M sh < 1O 1O M to compute 
their contribution to the total flux. 
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FIG. 1: Sky map of the annihilation signal from nearby 
subhalos of masses [10~ 6 , 10 6 Af Q ] contributing more than 
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FIG. 2: 3D view of the total $ c ° sm ° contributed by subhalos 
and the Galactic foreground, in a 100° x 100° f.o.v. centered 
on the GC. The Galactic foreground is cut at the subhalos' 
level. 



Figs. |21and|3]show the expected contribution to <j> cosmo 
accounting for both subhalos and Galactic foreground in 
the direction of the GC and the AC, respectively. As 
expected, the Galactic foreground is negligible in Fig. |3J 
while it is dominating around the GC in Fig. El Note 
that the very prominent peak due to the Galactic fore- 
ground at the GC position has been artificially truncated 
to appreciate the subhalos contribution. We estimate 
that $ c ° sm ° = 0.04 GeV 2 cm~ 6 kpc sr at the GC. 

Figs. [2]and|21also show that even the biggest contribu- 
tion to $ cosmo from a subhalo is well below the 0(1) level 
which, given our best value of $ SUSY ! would be required 
for detection in the GLAST experiment. In fact, in no 



FIG. 3: Same as Fig[2]in the direction of the AC. 



Monte Carlo realizations we have found that the subhalos 
contribution to $ cosmo exceed the value of O(10~ 3 ). 

To make this statement more quantitative and to 
also account for the possibility of detecting the 7-ray 
annihilation line we have evaluated the sensitivity of 
GLAST to the differential spectrum of photons expected 
from DM annihilation, convolved with its energy reso- 
lution AE—10%. We define the sensitivity a(AE) as 



n DM (AE)/ v / n hkg {AE) = 



T$e An 



Iae A f( E , 6 ) [d^ M /dEdn}dEdfl 



IAE 



bkj 



A off 



{E, 9)[d(j) bkg /dEdtt]dEdVL 



where T$ defines the effective observation time, 6aq = 0.7 
is the fraction of signal events within the optimal solid an- 
gle Af2 corresponding to the angular resolution of 0.1°, 
and A eS = 10 4 cm 2 is the effective detection area. We 
optimistically assume that the photon and charged par- 
ticles' detection efficiencies e 7 and e c h are 100%, thus 
we only consider galactic and extragalactic 7-ray back- 
ground as it is extrapolated by EGRET data at lower 
energies [28|. We simulated a continuous 5 years obser- 
vation of a pixel of our grid located at tp = 55°, assuming 
that in that pixel the value of <J> cosmo is the largest value 
found among all our Monte Carlo realizations. 

We chose a neutralino mass of 100 and 300 GeV, whose 
7-lines would in principle be observable with GLAST. 
The results are shown in FigQ] The black curves, corre- 
sponding to the left y-axis, show the sensitivity a(AE) 
of GLAST as a function of E, for two different values of 
the neutralino mass. The red curves, to be read on the 
right y-axis, show the expected flux from the same pixel. 
It is evident that GLAST would hardly detect neither 
continuum flux nor 7-lines from Galactic subhalos. 

Since it seems implausible to significantly increase 
$ SUSY , our results confirms those of OBElEl show 
that currently planned, satellite-borne experiments such 
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FIG. 4: Left y axis (black): Experimental sensitivity of 
GLAST to the photon flux from DM annihilation coming from 
a subhalo at ip ~ 55°, with $ c ° sm ° = 10" 3 GeV 2 cm" 6 kpcsr 
for m x = 100 GeV (solid line) and 300 GeV (dashed-line). 
Right y axis (red): Expected flux from the same subhalo as 
observed by GLAST with an energy resolution of 10%. 

as GLAST will not be able to detect the annihilation sig- 
nal produced by Galactic subhalos. This is due to the 
fact that, as pointed out by [15|. the total annihilation 
flux is dominated by massive Galactic subhalos rather 
than by minihalos. A result which makes our conclusion 
insensitive to variations in the low-mass density cut-off, 
in agreement with p^|. 

All plausible cosmological and astrophysical effects like 
the existence of a central core in the halo density profile, 



the various dynamical disturbances that reduce the sub- 
halo survival probability, the possible existence of a cutoff 
scale ~ 0.1fc cu t found [23 in the CDM power spectrum, 
generated by acoustic oscillations with wavelength com- 
parable to the size of the horizon at kinetic decoupling, 
would decrease the expected annihilation signal, making 
it our conclusions for subhalos more pessimistic. Unless, 
of course, one adopts steeper density profiles that, how- 
ever, are not supported by recent numerical experiments. 
Note that the 7-ray luminosity of our minihalos, which 
is designed to match the prediction of pj , is significantly 
smaller than that of j a difference that traces back to 
the large internal density of their minihalos. 

The presence of a minihalo population like that consid- 
ered in this work also enhance the 7-ray flux from nearby 
objects like M31. We estimate that, for a NFW profile, 
the total $ cosmo due to both the smooth and the subhalo 
contribution will not exceed the value of 5 x 10 -4 , making 
it impossible to detect. 

Recently, the HESS telescope has announced the 
serendipitous discovery of an unidentified extended TeV 
7-ray source, with a total flux above 380 GeV of ~ 
1CP 11 cm~ 2 s~ : 30]. If explained in terms of the annihi- 
lation of O(10 TeV) neutralinos, such a large flux could 
only be accounted for by advocating a central density 
r -1 - 5 , steeper than the NFW model. A prominent cen- 
tral spike might form around a supermassive black hole 
|3ll l32| that, according to hierarchical models of galaxy 
formation |33| , could reside at the center of both massive 
Galactic and sub-Galactic halos. 

We acknowledge fruitful discussions with L. 
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